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According to the Anderson theorem, the critical temperature T c of a disordered superconductor 
is determined by the average density of states and does not change at the localization threshold. 
This statement is valid under assumption of a self-averaging order parameter, which can be 
violated in the strong localization region. Stimulating by statements on the essential increase 
of T c near the Anderson transition, we carried out the systematic investigation of possible 
violations of self-averaging. Strong deviations from the Anderson theorem are possible due to 
resonances at the quasi-discrete levels, resulting in localization of the order parameter at the 
atomic scale. This effect is determined by the properties of individual impurities and has no 
direct relation to the Anderson transition. In particular, we do not see any reasons to say on 
"fractal superconductivity" near the localization threshold. 
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1. Introduction 

The general picture of coexistence of supercon- 
ductivity and the Anderson localizati- on was formed 
in the papers by Bulaevskii and Sadovskii [l]-[5] (see 
also [6j[7]). According to the Anderson theorem jS], 
the critical temperature T c of a disordered supercon- 
ductor is determined by the average density of states 
and does not depend on the form of one-particle 
eigenstates. Since the average density of states does 
not have singularity at the Anderson transition, so 
T c has the analogous behavior. The coefficient of 
the gradient term in the Ginzburg-Landau expan- 
sion, determining the superconducting response of 
the system, remains finite at the critical point. In 
the localized phase, the system breaks up into quasi- 
independent blocks of size £ (£ is the localization 
length) and superconductivity is suppressed due to 
the size effect, when the average level spacing in such 
a block becomes greater than T c . 

Recently it was stated by Feigelman et al (9[ [10] 
that T c increases at approaching the Anderson tran- 
sition from the metallic side and continues to grow 
in the localized phase: this is related with multi- 
fractality of wave functions. More than that, T c 
depends on the Cooper interaction constant g not 
exponentially, but in the power-law manner. For- 
mally, this statement does not contradict to the pa- 
pers [l]-[5]. Indeed, the Anderson theorem is valid 
under assumption of a self-averaging character of the 
order parameter, which in fact reduces to its spa- 
tial uniformity. According to estimates of |3J 0], 
the self-averaging property tends to violate when the 
localization threshold is approached and the space- 
inhomogeneous superconductivity is expected in the 
deep of the localization phase; so the true T c can be 
greater than its value given by the Anderson theo- 
rem. In fact, controversy between the papers [JJ-[5] 
and [9l [10] has an ideological character. The au- 
thors of p]-[5] proceed from the standpoint that lo- 
calization counteracts to superconductivity, so the 
latter encounters a lot of problems in the localized 
phase [2 [5] . Contrary, the results of [9[ [10] imply 
that superconductivity not only " survives" but even 
"prospers" in the localized phase. It looks suspicious 
from the physical viewpoint and contradicts the ex- 
perimental situation, which is in complete agreement 
with [TJ-[5]. 

The present paper has an aim to clarify a situa- 
tion. In fact, the essence of the problem is: how and 
in what extent self-averaging of the order parame- 



ter can be violated? The efficient approach to such 
problems was developed in [UJ-pT] and consists in 
the study of individual defects and their influence 
on the transition temperature. In particular, for the 
plane defects arranged perpendicularly to the z axis 
with the period L along it, the change of T c is de- 
termined by the formula^ 

ST 1 f 

Tf 1 = T37 9 2 \ dz [ v o v i 0) + v \ OO 2 ] . (!) 
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if there are no surface states localized near defects. 
Here g is the Cooper interaction constant, v\(z) is 
a deviation of the local density of states v(z) from 
its unperturbed value vq, \q = gv$ is the dimension- 
less coupling constant, T cu is the transition temper- 
ature in the absence of defects, integration is carried 
out over a vicinity of the single defect. For weak 
defects, only the linear in v\(z) term is essential, 
which exactly corresponds to the Anderson theorem 
and relates the change in T c with the change of the 
average density of states. Generally, v\ (z) is compa- 
rable with vq and already Eq.l predicts a possibility 
of essential violation of the Anderson theorem. It is 
related with the fact that the initially uniform order 
parameter is influenced by strong defects and can 
increase or decrease in their vicinity. More essential 
violations of the Anderson theorem are possible, if 
the surface states localized near defect appear at the 
Fermi level (Fig.l). In this case, the order parame- 
ter can be strongly localized near the plane defects, 
so T c does not depend on L and is determined by 
the BCS formula T c = 1.14o;o exp(-l/A 2 z)) with the 
coupling constant X2D, corresponding to the sepa- 
rated two-dimensional band (Fig.l). A crossover be- 
tween two regimes is appeared to be very sharp and 
the intermediate situation is of little interest. For- 
mally, the described results correspond to the pe- 
riodical arrangement of defects, but their character 
shows that the assumption on periodicity is not es- 
sential; so they give a complete picture for the small 
"impurity" concentration in the ID geometry. 

Analogous effects are possible in case of the point 
defects, where the localized regime for the order pa- 
rameter is related with existence of the quasi-local 
states (Fig. 2). A detailed investigation of these ef- 
fects allows to obtain the complete picture of possi- 
ble violations of self-averaging. The main conclusion 

1 We accept for simplicity that a plane defect changes only 
the density of states. Generalizations of (1), accounted for the 
change of the interaction constant g [121 1 13 1 and the cut-off 
frequency u>o |14| are also available. 
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Fi gure 1 : If plane defects are arranged perpendicular to the 
z axis, the problem allows a separation of variables and the 
longitudinal quasi-momentum fen can be introduced. For fixed 
fen , the spectrum is a set of bands with discrete levels splitted 
from them; if dependence on feii is taken into account, these 
levels turn into 2D bands, which can appear at the Fermi 
level. 

is that such violations are determined by individual 
defects and have no direct relation to the Anderson 
transition. One can distinguish two typical situa- 
tions. 

If disorder is created by weak impurities (Fig. 3, a) , 
then the assumption of self-averaging is always true 
and the Bulaevskii-Sadovskii picture is literally ap- 
plicable. The mobility edge lies near the initial band 
edge and T c is falling quickly at approaching it from 
the metal side due to decrease of the density of states; 
hence, superconductivity becomes practically unob- 
servable before the mobility edge is reached. Such 
situation is typical for the traditional superconduc- 
tors, which are good metals and effectively screen 
any impurity which is introduced in them. The ex- 
perimental situation is in complete agreement with 
these considerations [5]. 

In case of strong disorder (Fig.3,b), the mobil- 
ity edge can be located in the region of the prac- 
tically uniform density of statesd, so the T c value 
given by the Anderson theorem does not fall in ap- 
proaching the localization threshold. In fact, the 
true T c appears to be much larger and corresponds 

2 According to results by Zharekeshev I15| for the strongly 
disordered Anderson model, there is a wide plateau for the 
density of states in the center of band. 



Figure 2: A strong point defect inserted into site no of an 
ideal lattice leads to appearance of the local (-Ei) and quasi- 
local (Eo) levels. The latter corresponds to the maximum of 
the local density of states u(E, no). 

to localization of the order parameter at the small 
number of the "resonant" impurities, which produce 
the quasi-local states near the Fermi level. In ac- 
cordance with papers [9l [101 , T c depends on the in- 
teraction constant g in the power-law manner, but 
contrary to them, it has no essential dependence on 
the Fermi level position. It removes an illusion that 
localization "helps" superconductivity. In the vicin- 
ity of the true T c , observation of superconductivity 
is practically impossible due to a small fraction of 
the Meissner phase and negligible values of the crit- 
ical current. The easily observable quantity is the 
"bulk T c ", which can be defined theoretically as a 
transition temperature of the system with removed 
"resonant" impurities. Such T c is qualitatively (but 
not quantitatively) corresponds to the Anderson the- 
orem and the Bulaevskii-Sadovskii picture is con- 
firmed at this level. In such a case, T c obeys the 
"rectangular" dependence (Fig.3,b), which exponen- 
tially weakly deviates from the horizontal line near 
the mobility edge E c , and exponentially weakly de- 
viates from the vertical line near the endpoint E* of 
superconductivity. Such situation is typical for high 
T c superconductors, where coexistence of localiza- 
tion and superconductivity is easily observable [5]. 




* * * * 

fo> ' (b) 



Fi gure 3: (a) If disorder is created by weak impurities, self-averaging of the order parameter is not violated; T c strongly falls 
near the mobility edge due to decrease of the density of states, (b) In case of strong disorder, the true T c is determined by 
localization of the order parameter at rare resonant impurities, while the "bulk T c " obeys to the "rectangular" dependence. 



The estimate for the true T, 
T c ~ ga 



1 -< An J 



(2) 



(a is the lattice spacing, J is a bandwidth, and d 
is the space dimension) indicates that the "room" 
superconductivity is a widespread phenomenon. In 
fact, the growth of T c with the increase of Ao is 
bounded by a quantity wo/ 71 ", where luq is a cut-off 
frequency. For the phonon mechanism, such upper 
bound corresponds to the values already attained in 
high T c superconductors, and their further increase 
requires the use of higher frequency Bose excitations. 
In addition, the observation of true T c is probably 
possible only with the use of the scanning tunnel or 
squid microscopy |16) . 

According to the results of [9j [10] 



fJ 



l/ 7 



/ ~ (Tc/jy 



(3) 



where / is a portion of volume occupied by supercon- 
ductivity, and parameter 7 = 0.57 is related with a 
fractal dimensionality of wave functions. We do not 
deny the existence of the order parameter configura- 
tions, leading to results of type (3) (Sec. 3), but Eq.2 
corresponds to the higher value of T c . If supercon- 
ductivity is considered as a variational problem, then 



it is possible to say that our trial function is more 
successful than one in [TO]. The result (2) cor- 
responds to the rare peaks of the order parameter, 
which occur at the atomic scale and have a volume 
concentration ~ T c /J. Formally, these results cor- 
respond to Eq.3 with 7 = 1 and do not contain any 
information on multifractality; hence, there are no 
grounds to say on "fractal superconductivity" [10] 
near the localization threshold. 

2. Anderson theorem and inequalities for T c 

A basis for description of the spatially inhomo- 
geneous superconductivity is given by the Gor'kov 
equation for the order parameter A(r) 



A(r) = / K(r,r')A(r')d d r' 



(4) 



with the kernel K(r, r') in representation of exact 
one-particle eigenstates <p s (r) 



K(r,r')=gTj2J2 



<P* s (r)<Ps(r')>p* s , (r)yv(r') 
(e s - iu))(e s , +iu)) 



(5) 

where e s are eigenenergies (counted from the Fermi 
level), and summation occurs over the Matsubara 
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frequencies ui n = nT(2n + 1) with integer n. Fol- 
lowing de Gennes [17], we use the frequency cut-off 
| a; | < uj , which corresponds to the electron interac- 
tion 

V{ry;u) = -g6{u Q -\u\)8{r-r'), (6) 

which is strictly local and can be specified indepen- 
dently of one-particle eigenstates (in contrast to the 
momentum cut-off in the original BCS formulation, 
where interaction is defined by the matrix elements 
over plane waves). In the absence of magnetic ef- 
fects, eigenstates (p s {r) can be taken real and their 
orthogonality leads to the sum rule [T7] 

J K(r, r')d d r> = gv F {r) In , (7) 

where vp(r) = v(Q, r) is the local density of states 
v(e,r)=J2\Vs(r)\ 2 5(e-e s ) (8) 

at the Fermi level. 

The Anderson theorem follows from Eq.4 under 
assumption of a self-averaging order parameter, when 
A(r) and K (r, r') can be independently averaged 
over disorder. Since (A(r)} does not depend on r 
due to the spatial uniformity in average, the use of 
the sum rule (7) gives 

(A)= 9 ^)lnH^(A>, (9) 

and T c is given by the BCS formula, which contains 
the average density of states (vp). The latter does 
not change at the Anderson transition point, sug- 
gesting the analogous behavior for T c . More detailed 
information can be obtained, if Eq.4 is averaged over 
variable r 

(A) = g\n^^(Mr)A(r)) . (10) 

The function A(r) can be considered as positive^, 
and one has 

Vrmn (A) < (f F (r)A(r)) < v max (A) , (11) 

where v m in and v max are the minimal and maximal 
values of vp(r). It gives inequalities for T c 

l.Uuj Q e- 1/9U ™ n <T C < l.Uoj e- 1/9 "™"* , (12) 

3 For real ip s (r), the kernel K(r, r') is positive, since it 
can be written as flT^^ \Gu,(r, r')| (see Eq.5). The Cooper 
instability corresponds to the minimal characteristic value and 
the nodeless eigenfunction (the Entch theorem) 1181 . 



which can be also obtained from the known theo- 
rems of the matrix theory [13] [Sec. 2]. According to 
Eq.12, the power law dependence of T c on the cou- 
pling constant g [9] [TO] is impossible, if vp(r) has an 
upper bound v max . 

Near the Anderson transition, there are system- 
atic reasons for growth of the vp{r) fluctuations j3] 
S]. As noted in [5], the correlator (v(E+u, r) v(E, r')) 
at r = r' coincides with the Berezinskii-Gor'kov 
spectral density [TO] , which is determined by the dif- 
fusion pole with the observable diffusion coefficient 
D(uj,q) [20]: 

(v(E + u, r) u{E, r)> ~ Re f , ^ . 

J —iu) + D[u, q)q z 

(13) 

In the metallic phase, the static diffusion constant 
D(0,q) is real, so (vp(r) 2 ) diverges at the transition 
point as D -1 . In the dielectric phase, the analogous 
estimate can be obtained from the self-consistent 
theory of localization [21] by iteration of Eq.112 in 

D{u, q) = (-w)d(q) + Lo 2 di (q) , 

d(q)~e, diOd^M" 1 (14) 

(t is a distance to the critical point), so (vp(r) 2 ) ~ 
| | 1 and fluctuations grow symmetrically on two 
sides of the transition^. Estimations of the correla- 
tor (13) at the critical point based on multifractality 
of wave functions [10] suggest the dependence a; -7 
for u) — > 0; if divergency is cut off at the scale T c , 
then v max ~ vq{J /Tc) 1 / 2 and the maximum value 
T c - 5 2 /t allowed by Eq.12 is in a qualitative agree- 
ment with 110) . Consequently, if the upper bound 
for T c is realized in Eq.12, then it reaches the max- 
imum value at the transition point depending on g 
in the power law manner. 

However, the distribution of quantities y s (r) has 
the power law tails [10) and Eq.13 determines nei- 
ther the typical, nor the maximal value of vp{r). 
In fact, the given estimate for T c is not reached 
for weak disorder and is exceeded for strong dis- 
order. Formally, the approach of |10j is question- 
able due to replacement of matrix elements Mijki = 
J d d rip.i(r)ipj (r)(pk(r)(fi(r) by their mean values with 
averaging independently of the order parameter. 

More efficient approach is based on the study of 
effects from individual impurities, since it allows to 

4 According to the self-consistent theory, D ~ r in the 
metallic phase 1211 . 
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work with specific realizations of the random poten- 
tial and contains no problems of averaging. Intro- 
ducing one impurity after another, one can easily be 
convinced (Sec. 4), that unbounded values of vpir) 
can arise only from existence of quasi-local states 
(Fig. 2). The problem of quasi-local states has a gen- 
eral character. Indeed, one can imagine such fluctu- 
ation of the random potential, that a finite region of 
space is isolated from its environment by the high 
barrier; the corresponded discrete levels can have a 
very weak broadening and, appearing close to the 
Fermi level, can lead to unbounded values of vp{r). 
Such problems are discussed in the next section. 

3. Resonances at quasi-discrete levels 

Suppose that a system has a discrete spectrum 
and only one state is close to the Fermi level; then 
we can retain only one term in the sum over s, s' in 
(5): 

U u 

= 5 A(T)^(r)^(r'), (15) 
Then Eq.4 gives 

A(r) = X v 2 Q (r), 

X = gA(T) J ^ (r')A(r')d d r' . (16) 
and self-consistency of these equations determines 

l=gA(T)h, h = J vi{r)d d r. (17) 

Calculation of A(T) is possible without the cut-off 
frequency taken into account, since the sum con- 
verges at large ui: 

A(T) = rV - 1 =— tanh— . (18) 

For the exact resonance (eo = 0) we have A(T) — 
1/4T, so 

r c = 5/4/4, (19) 

and T c has a power law dependence on the interac- 
tion constant g. In the general case (see Fig. 4, a) 

T c = t—, : x °, , r , e c = 5/4/2 , 

ln(e c + e J - ln(e c - e ) 

(20) 



and a solution exists under condition 

|eol < e c . (21) 

At first glance, the considered regime is destroyed 
due to fluctuations!! or coupling with the continuous 
spectrum; in fact, it is not so (see below) and the 
main problem consists in the possibility to match 
the discrete level with the Fermi energy. 

Indeed, let the system has a finite size /, while 
its eigenstates are extended. Then the Fermi energy 
is located between two discrete levels^], and eo is 
determined by the average level spacing J(L/a)~ d ; 
estimating I4 ~ L~ d from the normalization condi- 
tion, we see that 

e ~J(L/a)- d , e c ~gL- d (22) 

and condition (21) cannot be fulfilled in the weak 
coupling regime, which is the only allowable in the 
BCS scheme. 

Let us couple the given system with a reservoir, 
and try to match the chemical potential of the lat- 
ter with the discrete level of the system. However, 
nothing good will occur from it: the local Fermi level 
of the system is still arranged between two discrete 
levels and it tends to equalize with the Fermi energy 
in the reservoir. The real flow of electrons is impos- 
sible due to elecroneutrality, and the problem will 
be solved by a minimal deformation: a double layer 
will arise between the reservoir and the system, and 
it will equate the Fermi levels. 

By the same reason, the situation cannot be im- 
proved due to localization of states. At first glance, 
in this case eo ~ J(L/a)~ d , e c ~ g£,~ d (£ is the 
localization radius of (fair)), so condition (21) re- 
duces to g > Ja d (L/£ t )~ d and can be fulfilled at suf- 
ficiently large L. In fact, blocks of size £ are quasi- 
independent and each of them has its own local Fermi 
level; these levels equalize due to double layers be- 
tween blocks, and the given estimates are valid only 
for L ~ £. In fact, the above arguments clarify the 
mechanism for the Coulomb gap [22] . 

It looks that the only possibility to avoid the 
given arguments is to take the size L of the atomic 
order. Indeed, at such a scale: (a) the notion of the 
Fermi level becomes senseless; (b) electroneutrality 

5 If eigenstate <p(r) is localized, then according to (20) a 
superconducting transition takes place in a finite system; of 
course, such conclusion is an artifact of the mean field theory 
and in fact the transition is destroyed by fluctuations. 

6 For a discussion of the parity effect see Footnote 16. 
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superconductivity 

Figure 4: (a) T c of the one-level system as a function of the level position eo in the absence of attenuation (solid line); 
attenuation 7 produces a shift of the curve by the quantity ~ 7 (dashed line), (b) The same, for a situation when the 
quasi-discrete level lies in the background of the continuous spectrum. 



can be violated; (c) a size of the double layer is com- 
parable with L. It means that the strong violations 
of the Anderson theorem can be exhaustively ana- 
lyzed by consideration of the one-impurity problem 
(Sec.4). 

Already at this stage it is possible to establish 
the relation with results of (9] [10] . In the considered 
there strictly one-electron picture, the discrete sys- 
tem of levels fluctuates freely relative to the Fermi 
energy, so resonances are possible at any length scale 
L. Then all principal statements of [HI [TO] are repro- 
duced: T c has a power law behavior as a function of 
g and does not depend on the cut-off frequency luq, 
while the order parameter A(r) follows the form of 
the wave function (see (16)) and will have multifrac- 
tal properties simultaneously with multifractality of 
the latter. However, this picture is completely de- 
stroyed, if electroneutrality is taken into account, 
since resonances at large scales become impossible. 
In fact, large scale fluctuations are insignificant even 
in a strictly one particle picture: a value of T c for an 
exact resonance, T c ~ gL~ d (see(17),(19)), is greater 
for small scales. 

Generally, the considered regime is not destroyed 
in the presence of the continuous spectrum. In this 
case, the level eg acquires the finite decay 7, which 



7 We have no doubt that papers [9l ll0| implicitly dealt with 
the same effect, but the improper averaging procedure led to 
a domination of large length scales. 



can be taken into account by replacement 

eo ± — > €q ±ioj ± 17 sign cj (23) 



A(T) = T 



1 



4 + iH + iY 



arctan ■ 

7re 7 + bT 



(24) 



where we have estimated the sum by the integral, 
introducing the cut-off |w| > bT (for the choice b = 
4/-7T such estimate practically coincides with the ex- 
act result (18) for 7 = 0). The finiteness of 7 leads 
qualitatively to the shift of the curve in Fig. 4, a by a 
quantity ~ 7, so a solution survives for 7 < e c . 
For finite Wq one obtains instead (24) 



A(T) » — arctan 

7T£0 



e (wo - bT) 



£2 + ( W0 +7)( 7 + &T) 



, (24') 



and can be easily convinced that finiteness of ujq is 
irrelevant under condition wo > e c . In the opposite 
case luq <C e c the allowed values of eo and 7 have 
an order (ecWo) 1 / 2 , while the maximal critical tem- 
perature T c is of the order ujq; in fact, restriction 
T c < ujq/tt is evident, since for T > luq/tt the sum 
over uj contains no terms. 

To investigate the effect of the continuous spec- 
trum on the order parameter, one can use the fol- 
lowing approximation for the kernel K(r, r') 



K(r, r') =K (r- r') + gA(T)ip 2 (r)<p 2 (r') , (25) 
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which ignores the backward influence of the discrete 
level on the continuous spectrum. According to [IT] 
I12j . such approximation provides qualitatively cor- 
rect description and can be justified in certain lim- 
iting cases. 

Having in mind a consideration of periodical con- 
figurations, we solve Eq.4 with the kernel (25) for a 
finite system of size L with the periodical bound- 
ary conditions. We accept L <C £o T ~ 1 ^ 2 I where 
r = (T — T c0 )/T c0 , £o is the coherence length, and 
T c o is a transition temperature, corresponding to the 
continuous spectrum. After the Fourier transform 
one has 



A q = gA(T)X 



1-K (q) 



(26) 



and self-consistency of two expressions leads to 



i = gA(T) 



iffl)q{vl)-q- 



-d 



Ko(q) 
-K (q) 



Using expansion in q 2 



1 - K Q {q) = A r + iAo£o¥ 



(27) 



(28) 



it is easy to see that the main contribution to the 
second sum in (27) occurs from small q, and the 
single term with q = is sufficient for L <C ^t^ 1 / 2 . 
Then Eq.27 accepts a form 



i = gA(T) 



h 



1 



(29) 



and the analogous approximations in (26) give 



A(r) = const 



1 



A i d 7 



(30) 



If T c o <C e C) then dependence of T c on eo has a 
form shown in Fig.4,b. In the zero approximation 
there are two independent systems, the quasi-local 
one with the transition temperature (20) (if attenu- 
ation 7 is small) and the continuous one characteriz- 
ing by T c o , while T c of the composed system is given 
by the maximal of two values. Interrelation of two 



8 In Sec. 4 we consider the one-impurity problem with the 
backward influence on the continuous spectrum. 



systems reduces to smoothing of dependence T c (eo) 
at the scale T c o(a/L) d / 2 , if ip(r) is localized at the 
atomic scale a. 

It is clear from (30) that the order parameter 
A(r) is practically constant for small r and local- 
ized at the scale a for large r. Crossover from one 
regime to another is very abrupt, and one can say on 
the "Anderson transition" for superconducting elec- 
trons. We see that the localized regime survives in 
the presence of the continuous spectrum, if the cor- 
responding T c exceeds T c q. In fact, existence of the 
continuous spectrum has a stabilizing effect on the 
localized superconductivity, since the order parame- 
ter takes non-zero values in the whole volume. 

4. One-impurity problem 

If G°„, is the Green function of an ideal lattice, 
V n = VS nno is an impurity potential, then the Green 
function G nn t of the perturbed system is determined 
by the Dyson equation [2"5] : 



G nn i — G nn i + G nng VG 



riQn' 



(31) 



Setting n — hq, one has the closed equation for 
G non > , whose solution is substituted into (31) 



G, 



Tin' ' nno rtnn' ' 



r 



nno ' nor 

V 



1 - VG° 



(32) 



where the scattering T-matrix reduces to a constant 
in the given case. For an ideal lattice G°„ does not 
depend on n, 



G° 



d d k 



(27r) d E - e k + iO 
is {e)de _ 



-I(E)-iTrv (E) (33) 

ti — e + iQ 

and so T-matrix has no no dependence. Condition 
1 — VI(E) = corresponds to existence of the local 
(if vq(E) — 0) or quasi- local (if vq(E) ^ 0) level 
(Fig.5). The local density at site n 



v(E,no) = 



ME) 



[1 - VI(E)] 2 + [kVv (E)}< 



(34) 



has an abrupt maximum near the resonance 1 
VI(eo) = (Fig. 6) with a value in it 



[v(E,n )} 



7r 2 ^o(eo) 



(35) 
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Figure 5: (a) If the impurity VS nno is inserted in an ideal lattice, equation 1 — VI(E) has two roots for large V, E\ and 
-Bo i which correspond to the local and quasi- local levels, (b) If the same impurity is inserted in the disordered lattice, both 
solutions correspond to the quasi-local levels. 
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Fi gure 6; The local density of states at the point no as a 
function of the impurity potential V. 



which grows unboundedly near the initial band edge 
(where fo(eo) — > 0). In the vicinity of the band edge, 
a calculation of G° n , is possible in the continual ap- 
proximation and gives at d = 3 (for \r — r'\ > a ): 

p ik F \r-r'\ 

G° nn , ^GV-r') = -™o-^-^, (36) 

Deviation of v{E, no) from v${E) is maximal for n = 
no and tends to zero for \n — n'\ — > oo. 

The Matsubara representation for the Green func- 
tions is obtained from (32) by replacement E — > 
ep + iu), where ep = for the corresponding choice 
of the energy origin. Setting no = 0, one has for the 
kernel in (4) 

^r,r')=.gT]T|G w (r,r')| 2 = 

UJ 

= 9Tj2K(r r')\ 2 + 

UJ 

+. 9 T]r Gl (r) % Gl{r') G°_Jr r')+ 

UJ 

+9TJ2 T- u G°_ u (r') Gl (r r')+ 

UJ 

+9TY, l G >)f l^l 2 K(rt = 
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K {r - r') + Ki{r,r') . (37) y u = f Ai (r) G°(r) d 3 r . (44) 



Solution of the Gor'kov equation with the kernel (37) 

is sought in the form Substituting Ai(r) from (436) into expressions (44), 



A(r) = A + A 1 (r) ! (38) 

where Ai(r) is localized near r — 0. Substituting in 
(4) and using the sum rule (7), one has 

A(r)= / K (r-r')A{r')d 3 r'+ 



and estimating arising integrals 

9 I Mr)\G°Jr)\ 2 d 3 r^X 01 



, . , 1.14wn . , , , 

+ 5i / 1 (r)ln^ r ^Ao + F(r) (39) 



|G°(r)| 2 |G«,(r)| 2 d 3 rEEA n , 
i/i (r) G° (r) d 3 r = A 02 + iA' ' 2 signw . 



where v\ (r) is deviation of the local density of states 

from vq and with the use of expressions for G° (r) and vi (r) (where 



/me real ana urn 
OAitrOdV. (40) and two primes 



.9 y G°(r)|G°,(r)| 2 d a r = Ai 2 +ai' 2S igna; (45) 

th the use of expressions for G° (r) and z/i (r) (where 
the real and imaginary parts are denoted by a prime 



Resolving (39) for A(r) by the Fourier transform and qO r r \ _ ni/ o CX p f _ M. r _|_ s jg n 

and simplifying the result analogously to (27), it is ^-F r I w f 

possible to separate the uniform term corresponding . 

to A , while the rest is attributed to Ai(r) ((. . ,) Vl ( r ) = I cos2fc f r + 7. sm2 ^ r _ ( 46 ) 



is the zero Fourier component): 



A ^0(°) 
A ° = TTT 



. , 1.14wn , . 

A ln^flWo + ^o 



A 1 (r) = A ln±^V 1 (r)+ F(r) . (41) 



(fcpr) 

it is easy to see that the integrals converge already 
for u, ui' = 0, so parameters Aoi, An, etc. can be 
considered as constant; it allows to write (44) in the 
form 

X^=X, Y U =Y' + iY"signu, 



Using the explicit expression for F(r) and setting in 

the integrals!! Z u = 2%Y' -2\T%\Y" + \Tjf X. (47) 



A 1 (r')G u (r-r')^A 1 (r')G° u (r), (42) 



The region remote from the resonance. The nat- 



r /• j i \ i ii r ural scale for the energy dependence of T-matrrx is 

one can transform (41) to the form . . °: f , 

given by the bandwidth J, so 7^ and |7^ | can be 

-Ko(O) 1.14wo considered as independent of ui anywhere, excepting 

A ° = A i 3 r ln T L9(^i)o A + ff(i/iAi) J , t he vicinity of the resonance (see below). Then Z u 

(43a) is also independent of us, and substitution of Ai(r) 

A / s A , 1.1 4wo / s n, 1^,(1/ M? from (436) into (44) leads to the linear system of 

Ai(r) = A ln — — — qvAr) + ql > 4 t" r . , . r . , „ , . .. , , , 

w T yiw equations for A and Z (see Appendix), whose sol- 

(436) ubility condition gives 

5T C _ 1 ^3 r ^i(r) + ^(r) 2 (4g) 



where 

= %Y U1 + T- u Y-u> + \%\ 2 X u , T c0 ~ A i 3 

f ,2 o Eq.48 is a natural generalization of the result (1): 

X = Ai (r) \C (r)\ d r 

u I iv J j\ ' the first term in the numerator corresponds to the 

Anderson theorem, while the second determines cor- 

9 This approximation is not quite rigorous but in fact it is rections to it . A configuration of the order param- 

used only tor estimates: the corresponding terms character- / \ 

izcd by parameters A 02 and A 12 have no significance both far eter shows that (48) corresponds to the delocalized 

from the resonance, and in its vicinity (see Appendix). regime. 
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For weak impurities (\V\ <C J) one has the esti- 
mates 



T' ra T ~ Ua 3 



T" ~ Ua 3 (U/J)(fc F a), (49) 



and 



^o(^i>0 ^ a 3 t/Q (V/J)(fef-a) 2 , 

(is 2 ) ~a 3 ^(V 2 /J 2 )(k F a)-\ (50) 

so the Anderson term is leading both in parameter 
V/J and in parameter (kpa)~ l . We accepted here 
fc F a <C 1, having in mind a situation near the band 
end, while estimates for the band center follow at 
k F a ~ 1. 

The delocalized regime retains in the case when 
the resonance condition 1 « VI{€f) is formally ful- 
filled, but the density of states v(e F ) is sufficiently 
large to provide a strong attenuation of the quasi- 
local state. In this situation 



T'-Ja 3 ^, T"~Ja 3 -. 

7 7 



(51) 



and one has under condition 7 > eo (where eo and 7 
are defined in Eq. 54) 



^0 (^i)o ~ (^1)0 ~ a v (k F a) 



(52) 



i.e. the Anderson term has the same order, as a 
correction to it. 

Vicinity of the resonance. If eo is a root of equa- 
tion 1 — VI(e), then in the vicinity of it 



1 - VI(e) = (e - e )/E 
and hence 

VE 



En ~ J . 



(53) 



T = 



7 = TrUT^M • (54) 



e - e + «7 
In the Matsubara representation one has 



iu; — eq + ij sign to 



Tj-iirj'lsignw, (55) 



where 



I'd = keo 



<3 + (M+7) a ' 
M+7 



^ + (l^l+7) 2 ' 



(56) 



so T-matrix can be considered as independent of u) 
under condition 



H<7. 



or 



(57) 



i.e. not very close to the initial band edge. If this 
condition is not fulfilled^, then the uj dependence is 
essential for the quantities 7^ and |7^'|, and hence 
for Z u . In fact, only one combination is relevant, 



(58) 



and substitution of Ai(r) from (436) into (44) al- 
lows to express through A and S; substituting 
this expression for Z^ into (58) and (43), one comes 
to the linear system of equations for Ao and S; its 
solubility condition with only leading terms retained 
(see Appendix) reduces to 



XnL 3 Vn 



-l + \ n V 2 E*A(T)) 



t z =Q. 



(59) 



where A(T) corresponds to expression (24). Equa- 
tion (59) describes the typical situation related with 
intersection of terms. The zero of the first bracket 
corresponds to the delocalized regime (see Eq.48), 
and vanishing of the second bracket corresponds to 
equation for T c of the localized superconductivity 
(compare with (17, 24)), while the term r 2 ~ (a/L) 3 
removes the degeneracy of terms in the intersection 
point (Fig. 7). 

5. Consequences for the Anderson model. 

Usually localization is studied in the framework 
of the Anderson model, which is a discrete version of 
the Schroedinger equation with a random potential: 
the bare spectrum is a band of width J, while the 
potential values V n on the lattice sites are indepen- 
dent random quantities with the distribution P{V} 
of width ~ W, which is supposed to be rectangu- 
lar. To transfer from the one-impurity problem to 
the Anderson model, it is sufficient to accept that a 
potential V of impurities fluctuates in the interval 
(— W, W), while their concentration c is gradually 
increased from small values to values of the order of 
unity. 

The results of Sec. 4 correspond formally to the 
periodic arrangement of impurities, but in fact their 
periodicity is not essential: each impurity arouses 
only a local deformation of the order parameter and 
these deformations are independent in case of a small 



10 In this case, the factor exp(— \u>\r/vp ) restricts contribu- 
tion to the integrals (45) by the atomic scale, where expres- 
sions (46) are inapplicable. 
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Figure 7: Dependence of T c on the impurity potential V for a small impurity concentration. 



concentration. If Ai(r) and fi(r) correspond to the 
one-impurity problem, then configurations 

A(r) =A + 5^A 1 (r-r i ), 

i 

i/p(r) = f + ^ v\{r - n) (60) 

correspond to a situation, when several impurities 
are arranged in points r% : it is a consequence of lo- 
calization of the kernel Ki(r,r') in both variables 
near the defect position. It is clear from (41) that 
the amplitude of Ai(r) is proportional to Ao, so 
Ai(r) = A /(r) and substitution of (60) into (10) 
gives for T c close to T c q: 

fH ^3 9jd 3 r[Mr) + Mr)f(r)} • (61) 

Here m is a number of impurities in the volume L 3 , 
and f(r) can be identified as fi(r)/vo from compar- 
ison with (48)0 One can sec that effect is propor- 
tional to a concentration of impurities, while their 
arrangement is irrelevant. The Anderson theorem is 

11 For weak disorder, relation Ai(r) = Ao^i(r)/^o follows 
from the second equation (41) after neglecting the quantity 
F(r), which is of the second order. Its validity for the delocal- 
ized regime without assumption of small v\ (r) is a non-trivial 
result expressed by equation (48). 



valid under condition |/(r)| <C 1, which is fulfilled 
for weak impurities. In case of nonequivalent impu- 
rities, the result (61) should be averaged according 
to distribution P(V). 

In the regime of the localized order parameter, 
each impurity is practically independent of enviro- 
ment and T c of the system is determined by those 
of them, which are close to a resonance; if distribu- 
tion P(V) is continuous and sufficiently wide, then 
the condition of almost exact resonance is always 
realized with a certain probability. Therefore, the 
concentration of the resonant impurities is finite and 
their quasiperiodic arrangement stabilizes the mean- 
field solution. 

Above considerations completely clarify a situ- 
ation for small impurity concentrations. Advance- 
ment to higher concentrations is simplified by ob- 
servation that equations (31), (32) never used a fact 
that G^ni corresponds to an ideal lattice; the same 
equations describe insertion of an additional impu- 
rity in the disordered superconductor. Noticing that 

nono ^E~e s + i0 J £ E-e + iO 

and replacing v(t, no) by its mean value (i/(e)), we 
obtain the same representation 1(E) — iitv^E) as in 
Eq.33, with a predictable behavior of 1(E) and v(E) 
(Fig.5,6). 
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Weak impurities. In this behavior of func- 

tions 1(E) and v(E) differs from their behavior in an 
ideal crystal by small smoothening of the Van Hove 
singularities (Fig. 5, b). Dependence on no results in 
fluctuations of the form of these functions, which 
are also small. It is clear that for weak impurities 
(| V\ -C J) the resonance condition is not fulfilled and 
no localization of the order parameter is possible. 

For the delocalized regime, it is convenient to 
present the result (61) in another form. Taking the 
one- impurity configuration A(r) = Ao + Aq f(r), 
v(r) = + Vx(r) and substituting it into equation 
(10), we have for the effective density of states en- 
tering into the BCS formula: 



y eff 



VQ + Mf) + ("l) + ("if) 

! + </> 



(62) 



Subtracting the result with / = and retaining the 
main terms in L~ d : 



y eff 



M = ("if) = 



(uif) v (ry 



AkpL 3 



(v) - vo = (v\) = - 



ty*T) 

2k 2 F L 3 ' 



(63) 



where we have taken into account that only the term 
with T' ~ V +V 2 I(Ep) is essential in Eq.46 for weak 
impurities. Inserting impurities one after another 
and averaging over V, 



ly eff 



(v) - vq 



W 2 

~ cvo — (k F a) 
W 2 

cvq -p (kFa) 



(64) 



we see that in the course of increasing a concentra- 
tion, the increment of the quantity v e ff — (v) is by 
a factor k F a smaller than the increment of (v) — uq. 
Near the band edge one have (v) ^> vq for sufficiently 
large concentrations, so v e ff ~ (v) ~ kpa(u) and de- 
viations from the Anderson theorem are small. Near 
the band center we have kpa ~ 1 and differencies 
(v) — vq and v e ff — (v) are small till concentrations 
c ~ 1; so v e ff — (v) <C (v). It is clear that violation 
of self-averaging for the order parameter does not 
occur for weak impurities. 

In the 3D case, isolated weak impurities do not 
produce bound states beyond the initial spectrum (it 
is clear from Fig. 5, a), and a finite density of states 
in this energy interval is a collective phenomenon re- 
lated with long-range fluctuations of the band edge. 



Consider a fluctuation in the region of size L, due to 
which the range of V values is somewhat restricted, 
(-W,W - 25) instead (-W,W). Then the mean 
value of the random potential is decreased by a quan- 
tity S, while a probability of such fluctuation exp(— L d 5/W) 
is not small for L d S/W < 1. Such fluctuations occur 
at all scales and produce a finite density of states 
beyond the bare spectrum. If size L of the fluctu- 
ation is sufficient for existing of superconductivity, 
the latter will not differ from superconductivity in 
the initial system with not shifted band edge; i.e. 
impurities will not violate the uniformity of the or- 
der parameter. Near the initial band edge, the in- 
dicated fluctuations strongly overlap and supercon- 
ductivity is quasi-homogeneous. Such fluctuations 
become spatially isolated in the region of strong lo- 
calization, where they can be described in terms of 
the size effect (Sec.6). 

Strong impurities. For a small concentration of 
strong impurities (\V\ > J), a behavior of functions 
1(E) and v(E) is not very different from their be- 
havior in an ideal crystal. However, in the region of 
the maximum of 1(E) the density of states becomes 
finite and, at first glance, complicates the occurrence 
of resonances. In fact, a new phenomenon comes to 
life. Since now G° oTlo depends on rig, the attenua- 
tion of the quasi-local state will be determined not 
by average density of states, but its local value at 
the point no, which can be small in a fluctuational 
manner. As a result, resonances become possible 
even for energies in the deep of the band, where they 
were forbidden in the ideal lattice. The typical situa- 
tion, when the local density of states ^f(^o) is small, 
corresponds to large values of the random potential 
in the vicinity of no; if now an impurity with large 
negative V is inserted into the site no, then a spe- 
cific resonance configuration arises (Fig. 8). El In the 
"minimal" variant, such configuration corresponds 
to existence of large barriers at the nearest neigh- 
bours of site no, while a value of the potential at no 
is chosen so that a corresponding level was in the 



12 The amplitude of long-range fluctuations can be seen 
from the fact that in the extremal cases the whole band is 
shifted by a quantity W or —W, i.e. such fluctuations by 
themselves (with no account for partial discretization of spec- 
trum) cannot produce unbounded values of vp(r). 

13 According to [25] , such configurations are responsible for 
multifractal statistics. It appears, that the tails of the dis- 
tribution function are determined by individual peaks (and 
not fractal clusters), in correspondence with our conception. 
Thereby, we do not ignore the existence of multifractality but 
give another description of its influence on superconductivity. 



13 




Figure 8: A typical fluctuation of the random potential re- 
sponsible for existence of the quasi-local state in the deep of 
the allowed band. 



interval of width ~ T c near the Fermi energy (the 
probability of this event is ~ T c /W). For a finite 
band, both large positive and large negative value of 
the potential are locking, and for W ^> J such values 
occur with probability p close to unity. Therefore, 
the probability of the "minimal" fluctuation 



(65) 



where Z is a number of the nearest neighbours. It is 
clear that such resonances can occur for any position 
of the Fermi level. In the region of the fluctuational 
tail, the density of states is small by the natural 
reasons and there is no need to create the barrier 
around no; so the factor p z will be absent but the 
less probable form of the effective potential well is 
necessary, in order the level was in the desired part 
of the spectrum. 

: 14 i 

With increasing of the impurity 
concentration, the effective bandwidth is extended 
and the maximum of 1(E) is shifted correspondingly. 
However, the general mechanism of resonances and 
estimation of their probability remain unchanged. 

Since the true critical temperature is hardly ob- 
servable, it is actual to consider the "bulk T c " , which 
can be defined as T c of the system with excluded res- 
onant impurities. For strong but not resonant impu- 
rities, two terms in Eq.48 are of the same order (see 
(50, 52)), and impurities are independent till concen- 
trations c~ 1, since the mobility edge lies far from 



14 Strictly speaking, the resonant configurations of such 
kind are possible for small W in the vicinity of the initial 
band edge. However, a size of such configurations is inevitably 
large (due to restriction of the barrier height and absence of 
levels in a shallow well of a small radius), so they have a negli- 
gible probability and are incompatible with electroneutrality 
(Sec.3). 



the bare edge of spectrum and kpa ~ 1. Validity of 
the Anderson theorem holds on the qualitative level: 
T c is determined by the effective density of states, 
which differs from the average one by a factor of the 
order of unity. 

6. Size effect in the localized phase. 

In the localized phase, the system breaks up into 
quasi-independent blocks of size £, and superconduc- 
tivity is suppressed due to the size effect. Below we 
analyze this effect in terms of the Gor'kov equation. 
Superconductivity in small samples was discussed in 
many papers (see a review article |26] ) , but this dis- 
cussion mainly concerns the aspects: 

(a) inadequacy of the grand canonical ensemble 
due to a fixed number of electrons in small granules; 

(b) parity effects; 

(c) insufficiency of the mean field approximation; 

(d) absence of an abrupt phase transition, etc. 
which are essential for finite systems and completely 
not actual in the present context. In principle, it is 
correct to stress unreliebility of the mean field ap- 
proach, but all attempts to overcome it (from modi- 
fied mean field approximations till the exact Richard- 
son solution and a direct numerical modelling) are 
based on the truncated BCS Hamiltonian, which by 
itself induces the certain way of pairing (in general 
incorrect)^- As for the Gor'kov equation, it cor- 
responds to the saddle-point approximation in the 
functional integral [27j [28] and is the most grounded 
of all mean-field type approaches; in addition, the 
electron interaction is specified in the physically clear 
manner and independently of one-electron states (Sec. 2). 
The accuracy of approximation is determined by the 
Ginzburg parameter, which provides insignificance 

of fluctuations in case of a superconductor (with ex- 
ception of some special cases: e.g. in finite systems 
fluctuations have a qualitative importance, destroy- 
ing a phase transition). The Gor'kov equation can 
be also obtained from the Eliashberg equations in 
the limit of the local interaction [M] . 

Consider the cubic sample of size L, accepting 
the periodical boundary conditions for the electron 
eigenfunctions. In a pure superconductor the latter 
have a form of plain waves, so |<p s (r)| 2 = L~ d and 
the local density of states (8) does not depend on 
r. Then A(r) = const is an exact solution of the 

15 The state ip s is coupled with its complex conjugated, 
which is correct only for a uniform order parameter I17| . 
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Gor'kov equation (4), which reduces to 



(66) 



and coincides with (9) in case of the continuous spec- 
trum. In a small energy interval, the spectrum can 
be considered as a set of equidistant levels with a 
spacing Q 

e a = n(s+l/2) 5 fl = l/v F L d , (67) 

where we accept that the Fermi energy lies in the 
middle between two discrete levels0. Substitution 
to (66) and summation over s gives 



1 = qT , F tanh 



O 



(68) 



For small f2, the argument of the hyperbolic tangent 
is large and one can set tanh a; = 1 — 2e 



-2.T 



J_ = ln l^o_ 4e - 2 ^T/n 
gv F T 



(69) 



where we retained only main terms with oj = ±ttT 
in the second sum over ui. Subtracting the analogous 
equation with 17 = 0, it is easy to obtain 



T c — T c q 



I _ 4 e -27r 2 T c0 /n 



(70) 



For T — > 0, one can replace summation in (68) by 
integration and obtain the equation for the critical 
value of f2, at which superconductivity is destroyed 



1 

gvF 
= In 



dbJ 1TLO 

— tanh — - 

2bJ ll r 



7TW0 



In a; 



■dx . 



(71) 



Jo cosh' 1 x 

The last integral is equal ln(7r/47), where In 7 = C = 
0.577 is the Euler constant and comparing with the 
result for T c q 



1 = gvp In 



47W0 



. , 27CJ0 , s 

1 = gu F In , (72) 



one can sec that 



fir 



2ttT, 



o0 • 



(73) 



16 Such assumption is commonly accepted |26| for the case 
of the even number of electrons N; for odd N it is accepted 
e s = Qs, but the level eo = is considered as "blocked", 
i.e. occupied by the unpaired electron and not participating 
in the scattering process. In the latter case, the results are 
analogous but correspond to smaller T c . 



To find the dependence of T c on Q, in the vicinity of 
f2 c , one can transfer (68) using the Poisson summa- 
tion formula [2"9"1 



1 00 



TTOJq/Q, 



gv F 



tanh a: f s£l 

dx exp < i— -x 

2x \ nT 



-TTLUa I f2 



(74) 

where the term with s = corresponds to (71). For 
s / 0, the integrals are convergent at large \x\ and 
it is possible to set ujq — 00 in them. Due to even- 
ness in s they can be calculated for s > 0; then the 
contour is shifted in the upper half-plain and the 
main contribution arises from the pole x — in/2. 
For Q/T ^> 1 it is sufficient to retain the terms with 
s = 0,±1, 



J_ = l n ±^_ 4e -^. 
gvF il 



(75) 



and subtracting the analogous equation with T = 



T = 

c 21n[4f2 c /(fi c -n)] ' 

In the reduced coordinates 



n -> n c (76) 



V = T c /T c 



c0 : 



x = o/ri c 



(77) 



one can obtain the universal dependence y(x). In- 
deed, transforming (68) by subtraction of the anal- 
ogous equation with Q = 0, one has 



In — = T V ^- (tanh 



7r|w| 



1 



(78) 



where wo can be set to infinity. Substituting the 
Matsubara values 7rT(2n+ 1) for lj, one can present 
the dependence y(x) in the parametric form 



y = exp F(t) , 



t~ 1 expF(t) 



1 



n=0 



2n + l 



tanh 



w(2n + l)t 



1 



(79) 



where t runs from zero to infinity. Numerical calcu- 
lation based on (79) gives the "rectangular" depen- 
dence y(x) shown in Fig. 9: this dependence has ex- 
ponentially small deviation from the horizontal line 
near y = 1, and exponentially small deviation from 
the vertical line near x = 1. 
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Figure 9 : " Rectangular" dependence of T c on the level spac- 
ing Q in a finite system; it is universal in the reduced coordi- 
nates y = T c /T c o, x = Q/Q c . 



The given consideration retains for a disordered 
superconductor if possibility of self-averaging is ac- 
cepted. The obtained results can be used to de- 
scribe the dependence of T c on the distance to the 
mobility edge in the localized phase, where the sys- 
tem is divided into quasi-independent blocks of size 
£. The role of f2 is played by the quantity 



n(E)~ j(e/«r d ~ j{\E-E c \/jy v . 



(80) 



where v is the critical exponent of the localization 
length. According to Sec. 5, the assumption of self- 
averaging is valid literally for weak disorder and on 
the qualitative level for strong disorder in the ab- 
sence of resonances. In the latter case, T c is deter- 
mined by the effective density of states which dif- 
fers from the average one by a factor of the order 
of unity, which is a smooth function of parameters. 
It preserves the character of singularities (70) and 
(76), which determine the behavior near E c and E* 
(Fig. 3) and are responsible for the most striking fea- 
tures in the dependence T c (cf)- 

7. Conclusion. 

The present paper resolves contradiction between 
two series of papers [T]-[5] and [5JED]- The obtained 
results has in some way a compromise character. On 
the one hand, the " bulk" superconductivity behaves 
in correspondence with the picture by Bulaevskii and 



17 Of course, in this case one should take some realistic 
statistics of the Wigner-Dyson kind instead of the equidis- 
tant levels, but it has a small effect on the results I26| . 



Sadovskii Q]-[5]. On the other hand, the true tran- 
sition temperature T c of strongly disordered super- 
conductor does not coincide with the "bulk" one and 
is determined by rare peaks of the order parameter 
on the atomic scale; in correspondence with [5J [TU] 
it has a power law dependence on the coupling con- 
stant and does not depend on the cut-off frequency. 
However, in contrast to 0110]) it has no essential de- 
pendence on the position of the Fermi level and does 
not correlate with the Anderson transition. By this 
reason, we do not see any grounds to say on "fractal 
superconductivity" [10] near the localization thresh- 
old. 

The obtained results are valid in the framework 
of the mean field theory and can be partially mod- 
ified by fluctuations. It is essential, however, that 
the role of fluctuations is determined by specific val- 
ues of parameters, characterizing the system: if, for 
example, the ratio T c /J is not too small, then the 
resonant impurities have rather large concentration 
and the Josephson coupling between the localized 
superconducting " drops" is sufficiently large for sta- 
bilization of the mean-field solution. Analogous re- 
marks are valid in relation with the Coulomb blocade 
effects [27]. 

The Anderson model with several types of peri- 
odically arranged inpurities can be considered as a 
model for high T c superconductors. A possibility to 
interprete the "superconducting explosion" of 1987 
as localization of the order parameter was indicated 
in [TT] ; the above results show that such localization 
is possible not only at the Cu — O planes but also 
at the individual atoms. The adequacy of such a 
model is confirmed by the following: (a) the maxi- 
mal T c values practically coincide with uiq/tt; (b) the 
estimates for T c are quite optimistic; (c) the isotope 
effect is practically absent in the regime e c <C loq- 

Appedix. On solution of the Gor'kov equation with 
the kernel (37). 

Let fill in the gaps for our exposition in the main 
text. 

In the region remote from the resonance, we can 
consider 7^ and |7^'| as independent of lo: then 
is also constant. Substituting Ai(r) from (436) into 
expressions (44) for X u and Y u , we have representa- 
tion (47) with parameters 

1.14wn UJn 
X = A \n — — Aoi + —XnZ 

1 7T 



Y' = A In 



Tf, A 02 H A 12 Z 

1 TT 
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r =A ln A 02 + — A 12 Z. 

Then Z,,, has a form 



1 7T 



(Al) 



(A2) 



and its combination with (43a) gives a system of 
equations for Aq and Z 



A 



1 I \ 2/ 2\ i 1-14^0 



+ 5 — A O iZ = 

7T 



A B 2 In 



1.14w 

r 



z 



U) 



B, 



0, (A3) 



with the coefficients 



„ Jfo(0) . 1.14w 
£>i = - — in ■ 



T 



B 2 — 2T'A 02 — 27 - "A ' 2 



|T| 2 A 



B3 = 2r'Ai 2 -2|r"|A' 1 ' 2 + |r| 2 A 11 . 



(A4) 



The terms containing ujq has an additional smallness 
~ ujq/J and can be neglected^; the condition of 
solubility for (A 3) gives the result (48). 

In the vicinity of the resonance, one cannot ne- 
glect the oj dependence of the quantities 7^, |7^'|, 
and consequently Z u . Substituting Ai(r) from (436) 
into (44) for and Y u , one has representation (47) 
with parameters 

X = Aolni^Aoi + AnS 



r' = A lni^A 02 + A' 12 5 



y =A ln - A 02 +A 12 5 



(A5) 



and for Z^ 



Zuj = A In l^p- [2%X', 2 - 2\TJ\ K 2 + 

-\%\ 2 \oi] + S [2%\' l2 - 2\TM + \%\ 2 Xn} 

(A6) 



18 We have in mind the traditional superconductors. If (Jo ~ 
J, then the "vicinity of the resonance" is extended and in fact 
occupies the whole band. 



Substitution into expressions (58) and (43) gives a 
system of equations for Aq and S 



A 



-^-+5(^1)0+5 (^i)ohi — — 

+gKiS = o, 
Aodhx^p^+S [-1 + C 2 ] = 0, 
with dehnitions 

G\ = 2A 02 (Ti — 2A ' 2 (T 2 + A01CT3 , 
C 2 = 2A' 12 (Ti - 2A 12 o"2 + A11CT3 , 



(A7) 



(A8) 



^3 = T^ \%\ 2 . 
The condition of solubility for the system (A. 7) gives 
(^-5(-i)o-3 2 K 2 )oln^) (-1 + C 2 ) + 

1.14wn 

+ C l5 A i In — ^ = (A9) 
Estimations for eo ~ 7 ~ T give 

(Ti - K-Efl , CT 2 ~ VEq In — , CT 3 ~ , 

(A10) 

and allow to retain only the leading terms in J/T; 
as a result, Eq.A9 can be written as 



Agl 5 ' (- 1 + Xll(T 3) + 



\2 

A^ 



5 CT 3 = 



(All) 



and reduces to a form (59); the last term is essential 
only near the intersection point of dashed lines in 
Fig. 7, when eo ~ 7 ~ e c and T should be replaced 
by e c in the estimates (A 10). 
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